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ABSTRACT 



We have analysed an XMM-Newton observation of the low mass X-ray binary and atoll source MXB 1728-34. The source was in a 
low luminosity state during the XMM-Newton observation, corresponding to a bolometric X-ray luminosity of 5 x 10 36 , kpc erg s . 
The 1-1 1 keV X-ray spectrum of the source, obtained combining data from all the five instruments on-board XMM-Newton, is well 
fitted by a Comptonized continuum. Evident residuals are present at 6 - 7 keV which are ascribed to the presence of a broad iron 
emission line. This feature can be equally well fitted by a relativistically smeared line or by a self-consistent, relativistically smeared, 
reflection model. Under the hypothesis that the iron line is produced by reflection from the inner accretion disk, we can infer important 
information on the physical parameters of the system, such as the inner disk radius, R m = 25 - 100 km, and the inclination of the 
system, 44° < i < 60° . 
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1. Introduction 



■ line: identification — stars: neutron — stars: individual: MXB 1728-34 — X-ray: binaries — X-ray: 



Broad iron emission lines in the energy range 6.4-6.97 keV 
have been detected in high energy resolution spectra of many 
X-ray sources c ontaining a compact object, such as act ive galac- 
tic nuclei (e.g., iTanaka et al.lll995t iFabian et al J 120001) and X- 
ray binary system s containing a stellar-mass black hole (e.g., 
iMiller et all 120021; iMillerl 120071. for a review), or a weakly 
magn etized neutron star (e.g.. iBhattacharvva & Sttohm avei 
2007; ICackett et ail 120081: iDi Salvo et alj 120091: iPapitto et al 
l2009tllaria et alj2009HD'Ai et al.l2009l and references therein). 
Identified with fluorescent Key transition of iron at different ion- 
ization states, these lines are generally interpreted in terms of 
refle ction of the central hard X-ray emission on the accretion 
disk (Fabi an et al.lll989 l). Under this hypothesis, these lines are 
made broad and asymmetric by Doppler and relativistic effects 
induced by the Keplerian motion in the accretion disk near the 
compact object. The shape of the line is therefore an almost 
unique proxy o f the innermost accretion d isk close to the com- 
pact object (see [Reynolds & Nowakll2003l for a review), and, in 
particular, on the inner disk radius. It also indicates the inclina- 
tion angle of the system and the ionization state of the reflect- 
ing matter. Other reflection components like absorption edges 
and the Compton hump, usually observed between 20-40 keV, 
are also expected resulting from photoelectric absorption and 
Compton scattering of the main Comptonization continuum on 
the accretion disk matter. 

MXB 1728-34 (4U 1728-34, GX 354-0) is a low mass 
X-ray binary containing a weakly magnetized accreting neu- 
tron star. The optical counterpart of this "galactic bulge" 
source has not been identified yet, due to the high optical 



extinction towards the Galactic center. Di scovered in 1976 
with the Small Astro nomy Satellite SAS-3 dLewinet all 1 19761: 
Hoffman et al. 1976 ), this source bel ongs to the so-called atoll 
class dHasinger & van der Klisll 19891) and shows frequent type- 
I X-ray bursts (e.g., Basinskaet al. 1984) which are caused by 
thermonuclear flashes on the neutron star surfac e. Furthermore, 
doubl e-peaked burst profiles have been observed (Hoffma n et al.l 
1976) which are explained a s caused by photospheric radius 
expansion during the burst (Taam| Il982l) . These bursts have 
been used to constrain the distanc e to the source, betwee n 4. 1 
and 5.1 kpc (IDi Salvo et alJ I200T3: iGailowav et all l200l . The 
power spectrum of MXB 1728-34 displays kilohertz quasi- 
periodic oscillations (QPOs) in the persistent emission, and a 
nearly coherent oscillation at ~ 363 Hz during some bursts that 
has been interpreted as the spin frequency of the neutron star 
(IStrohmaver et al.lll996b . 

Spectral analysis of MXB 1728-34 has been performed in 
the past using dat a from different sat ellites, such as Einstei n 



EXOSAT (Whi te etal 
ROSAT dSchulz! ' 



(iGrindlav & Hertzl [l98ll). SAS-3 dBasinska et all 11984), 
19861) . SIGMA dClaret et all 1 19941) . 
and more recen tly using BeppoSAX 



2000), ASC A dNarita et al.l 
200f4_ INTEGRAL 



1999) 



(Pirai no etal.1 120001: IDi Salvo et alj 
120011). 



RXTE, Chandra (D'Aietal 



mail to elise.egron@dsf.unica.it 



dFalanga et aT]l2006l) . and XMM-Newton dNg et al. 2010b- The 
X-ray spectrum is generally composed of a soft and a hard 
component. The first one can be described by a blackbody or a 
multicolor disk blackbody, and may originate from the accretion 
disk. The second one can be fitted either by a Comptonized 
spectrum or a thermal Bremsstrahlung. The Comptonized model 
seems to be more realistic, since it is efficiently produced by soft 
photons coming from the neutron star surface and/or boundary 
layer between the accretion disk and the neutron star and/or 
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Table 1. Instrument modes, filters and exposure times. 



Instrument 


Mode 


Filter 


Exposure (ks) 


pn 


Timing mode 


Thick 


26.9 


MOS1 


Timing mode 


Thick 


27.5 


MOS2 


Timing mode 


Thick 


27.5 


RGS1 


Standard spectroscopy 




28.1 


RGS2 


Standard spectroscopy 




28.1 



the accretion disk, which are up-scattered by electrons in a hot 
corona. 

A broad emission line at 6.7 keV has been often detected 
in the X-ray spectra of this source and h as been interpreted a s 
emission from highly ionized iron (e.g., iDi Salvo etaf] 120001) . 
The large width of the line s uggests that it cou ld come from an 
ionized inner accretion disk (Pirain oet al.l20 00). or alt ernatively 
it coul d be emitted from a strongly ionized corona. iD'Ai et alj 
(2006) have proposed an alternative model to describe the iron 
line region using two absorption edges associated with ionized 
iron instead of a Gaussian line. 

In this paper, we present a spectral analysis of high energy 
resolution data taken by XMM-Newton on 2002 October 3rd us- 
ing all the five X-ray instruments on-board t his satellite. Th e 
EPIC-pn data have already been published by iNg et al.1 (1201 Oh . 
in a 'catalog' paper dedicated to the study of the iron line in 
16 neutron star LMXBs observed by XMM-Newton. Here we 
present a different approach to the analysis of these data, since 
for the first time we fitted XMM-Newton data of MXB 1728-34 
with a self-consistent modelling of the continuum emission and 
of the reflection component. We also tried several models to fit 
the iron line profile, such as diskli ne or relline, wh ich are differ- 
ent from the models proposed by Ng et al. (2010) (Gaussian or 
Laor). All our results favor the relativistic nature of the line pro- 
file. We fitted simultaneously the spectra from all the five X-ray 
instruments onboard XMM- Newton, while just the pn data are 
analysed by INg et al.1 (1201 Oh . 



2. Observation and data reduction 

MXB 1728-34 was observed by XMM-Newton on 2002 October 
3rd for a total on-source observing time of 28 ks. The ob- 
servation details for the instruments on-board XMM-Newton 
including the Eur opean Photon Imaging Camera (EPIC-pn, 
Struderetal.120011) , the MOS 1 and MOS2 cameras dTurner et al.1 
2001), and the Reflection G rating Spectrometer (RGS1 and 



RGS2, Iden Herder et al .2001) are presented in the Table 1 . The 
Optical Monitor (OM, Maso n et al JI200 lb was not active during 
this observation. 

The 1.5-12 keV lightcurve of the All Sky Monitor (ASM) 
on board RXTE extracted within ~ 450 days the XMM-Newton 
observation indicates that the source was not in a high activity 
state, since it shows an average count rate of about 3 counts/s. 

The XMM-Newton data were processed using the Science 
Analysis Software v. 9. The EPIC-pn camera was operated in 
timing mode to prevent photon pile-up. We created a calibrated 
photon event file using the pn processing tool epproc. Before 
extracting the spectra, we checked for contamination from back- 
ground solar flares by producing a lightcurve in the energy range 
10-12 keV. There were no high background periods during this 
observation. We used the task epfast to correct rate-dependent 
CTI effects in the event list. The source spectrum was extracted 
from a rectangular area, covering all the pixels in the Y direction, 
and centered on the brightest RAWX column (RAWX=38), with 




52300 52350 52400 52450 52500 52550 52600 52650 52700 52750 
Time (MJD) 

Fig.l. RXTE/ASM lightcurve covering the 1.5-12 keV energy 
range. The vertical line indicates the time of the XMM-Newton 
observation performed on 2002 October 3rd. 



a width of 13 pixels around the source position (because 90% of 
the source counts up to 9 keV is encircled by 53 arcsec, and 1 pn 
pixel is equivalent to 4.1 arcsec). We selected only events with 
PATTERN < 4 (single and double pixel events) and FLAG=0 as 
a standard procedure to eliminate spurious events. We extracted 
the background away from the source (in the RAWX = 6-18). 
We also checked that pile-up did not affect the pn spectrum us- 
ing the task epatplot. The total count rate registered by EPIC-pn 
CCDs was around 1 10 count/s, and was 64 count/s in the 2.4-1 1 
keV range, slightly increasing (by a 5%) during the observation. 

The MOS data were also taken in timing mode and processed 
with the routine emproc to produce calibrated event list files. The 
source spectra were extracted from a rectangular box centered on 
RAWX=320 (MOS1), and on RAWX=308 (MOS2), selecting 
an area 30 pixels wide around the source position, and covering 
722 pixels on the Y (PHA) direction. Only events corresponding 
to PATTERN < 12 and FLAG= were selected, correspond- 
ing to standard filters. The background spectra were extracted 
far from the source, centered in the column RAWX= 240. We 
checked that the MOS spectra were not affected by pile-up. The 
count rates were estimated around 30 counts/s for each MOS 
unit (20 counts/s considering 2.4-1 1 keV energy range). 

Spectral channels of EPIC-pn and MOS spectra have been 
rebinned to have 3 channels per energy resolution element and 
at least 25 counts per energy channel. 

The two RGS were operated in the standard spectroscopy 
mode. RGS data were processed using the rgsproc pipeline to 
produce calibrated event list files, spectra and response matrices. 
The count rates measured by RGS1 and RGS2 were around 2.5 
and 3.5 counts/s, respectively. The RGS data were rebinned in 
order to have a minimum of 25 counts per energy channel. 

3. Spectral analysis 

Data were fitted by using Xspec (lArnaudl 19961) v. 12.5.1. All un- 
certainties are given at the 90% confidence level (Ay 2 = 2.706). 
We fitted simultaneously the broad band energy spectra of the 
source obtained from all the five instruments. Considering the 
best calibration ranges of the different detectors, the data analy- 
sis from EPIC-pn, MOS 1 and MOS2 cameras was restricted to 
the energy range 2.4-11 keV. This excluded the region around 
the detector Si K-edge (1.8 keV) and the mirror Au M-edge (2.3 
keV) that could affect our analysis. This problem was already no- 
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Table 2. Best fitting parameters of the continuum emission for 
the XMM-Newton pn, MOS1, MOS2, RGS1 and RGS2 spectra 
of MXB 1728-34. 



Parameter 



Value 



Nh (xlO I2 cm- z ) 
£T seed (keV) 
kT e (keV) 

T 

Norm 

Flux 2.0-10.0 keV (pn) 
Flux 2.0-10.0 keV (MOS) 
Flux 1.0-2.0 keV (RGS) 
Total x 2 (d.o.f.) 



2.4 ±0.1 
0.59 ± 0.02 
2.74 ± 0.04 
16.5 + 0.2 
9.54 ± 0.02 
8.06 

8.17 (MOS1) - 8.07 (MOS2) 
0.186 (RGS1)- 0.185 (RGS2) 
1732 (903) 



Note: The model used to fit the continuum is cons*phabs*compTT. The 
absorbed flux is in units of 10~'° erg cmT 1 s _1 . 



ticed for the EPIC-pn o b servations performed in timing mode , 



Boirinetal. (2005), Iariaetal. (2009), D'Aietal 



Papit to et alJd2009l) . lD , Ai et alJ(l2010l) . lPapitto et alJd2010h . We 



used only the RGS1 and RGS2 data between 1-2 keV to con- 
strain the softest band. 

The different cross calibrations of the five instruments were 
taken into account by including normalizing factors in the 
model. These factors were fixed to 1 for pn and kept free for 
the other instruments. 

We first fit the cont inuum wi t h a th ermal Comptonized 
model using compTT (iTitarchukl 1 19941) . modified at low 
energy by the interstellar photoelectric absorption mod- 
elled by phabs using photoe l ectric cross-sections of 
Balucinska-Church & McCammon (1992) with a new He 
cross-section based on Yan et al.l (1 19981) and standard abun- 



dances of lAnders & Grevessd (Il989l) . The ^/degrees of 
freedom (d.o.f.) of the fit was large, 1732/903. We then tried to 
add a blackbody component (bbody model) to improve the fit. 
The addition of this component resulted to be not statistically 
significant, thus we decided not to include the blackbody in our 
model. The values of the parameters of the continuum emission 
are reported in Table 2. 

With respect to this continuum model, an excess was present 
in the residuals between 5.5 and 8 keV, probably indicating 
the presence of iron discrete features. The fit was improved 
by adding a broad iron emission line, modelled by a simple 
Gaussian line (Model 1 in Table 3), centered at 6.6 keV with the 
cr parameter frozen at 0.6 keV. Indeed, the width of the line was 
poorly constrained, however when we fitted the MOS data inde- 
pendently from the pn data, a significantly broadened line was 
detected, with a lower limit on cr of 0.55 keV. With the addition 
of the Gaussian, the fit gave a^ 2 /d.o.f. = 1489/901 (resulting to 
a significant improvement of the fit, with a Ax 2 = 243 for the 
addition of two parameters). We also tried to fit the line with a 
combination of two emission lines instead of a broad line, but the 
X 2 was worse (x 2 /d.o.f. = 1531/901). So this broad line does not 
result from a blending of iron line at different ionization state. 

We then tried to substitute the Gaussian at 6.6 keV with 
a diskline profile and obtained a slightly better result with a 
^ 2 /d.o.f. = 1463/899 . The data and the residuals obtained us- 
ing this model are shown in Fig. 2. Due to the large uncertainties 
on the outer radius of the disk and on the inclination of the sys- 
tem if they were let free, we have frozen the first one at 1000 
R g and the inclination at 60° (the source does not show any dip 
in its lightcurve, implying i < 60°). The improvement of the fit 
corresponds to Ax 2 = 26 for the addition of two parameters (the 
F-test gives a probability of chance improvement of about 10~ 4 ). 



We think that th e F-test gives a reliabl e result in this case. 
In fact, according to Protass ov et al.l (120021) . the conditions that 
have to be satisfied in order to properly use the F-test are: i) 
the two models that are being compared must be nested; ii) 
the null values of the additional parameters should not be on 
the boundary of the set of possible parameter values. We think 
that a comparison between a Gaussian and a diskline satisfies 
these conditions, since a Gaussian profile can be obtained by a 
diskline (which is indeed given by a narrow, cr = 0, Gaussian 
multiplied by the kernel rdblur) and a good approximation of a 
Gaussian profile is obtained by a diskline for values of the pa- 
rameters not at their boundary. However, in order to assess the 
significance of the relativistic line smearing on a statistical ba- 
sis, we used an other statistical m ethod based on the Posterior 
Predictive p Values described by Hur kett et al.l J2008). For sim- 
plicity, we have restricted our data to two instruments (pn and 
MOS2). These data were fitted with Model 1 and Model 2, re- 
spectively. We obtained an improvement of the fit correspond- 
ing to a Ax 2 of 18 for the addition of 2 parameters when we 
substituted Model 1 with Model 2. Then, we simulated 200 pn 
and MOS2 spectra according to Model 1, which were fitted in 
a second step using Model 2. The Ay 2 was registered for each 
simulation. Among the 200 simulations, we found 2 times a Ay 2 
higher than 18. So the probability of chance improvement we 
get from these simulations is 1%, which is in agreement with 
the 0.75% calculated by the F-test performed on the real data re- 
stricted to two instruments. Therefore we can conclude that the 
diskline model is preferred since it gives a probability of chance 
improvement of the fit of about 10~ 4 using all the instruments. 

This model (Model 2 in Table 3) gives an estimate of the 
inner radius of the disk R m ~ 18 R g (R g = GM/c 2 is the 
Gravitational radius). Fixing the inclination to lower values (i.e. 
i < 60°) , worse x 2 were obtained with the other diskline pa- 
rameters drifting towards lower R m , higher rest-frame energies, 
and higher, in absolute value, emissivity indexes. We searched 
for an absorption edge in the energy range 7-10 keV, but none 
was significantly detected. 

To fit the iron line, we also used a new model for a relativis - 
tically distorted disk line, called rellin^H (Dau ser et alj 12010). 
which calculates line profiles taking into account all the rela- 
tivistic distortions in a disk around the compact object. We also 
fixed the outer radius of the disk and the inclination of the sys- 
tem to the same values used in the diskline model. The^- 2 /d.o.f. 
obtained in these conditions is 1464/899. The best fit line param- 
eters obtained in this way are perfectly consistent with those ob- 
tained using the diskline model. In particular this model (Model 
3 in Table 3) estimates the inner radius of the disk to be Ri n ~ 19 
R s . 

We tried an alternative model for the iron features (see 
iD'Ai et aTl 120061) using two absorption edges (instead of an 
emission line) which are found at 7.50 keV (t ~ 0.06) and 8.49 
keV (t ~ 0.06), associated to mildly and highly ionized iron, re- 
spectively. The^ 2 /d.o.f. for this fit is 1519/899 (which has to be 
compared to 1489/901 that we obtained fitting the iron line with 
a Gaussian (Model 1) or to 1463/899 that we obtained fitting the 
iron line with a diskline (Model 2)). Therefore this model (called 
Model 4 in Table 3) gives a worse fit of the iron features than the 
previous ones. 

In order to test the consistency of the broad iron line with a 
reflection component, we fitted the data using reflion, a self- 
consistent reflection model including both the refl ection con- 
tinuum and the corresponding discrete features (Ross & Fabian 
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Table 3. Results of the fitting of MXB 1728-34 XMM-Newton data with different models. The five models reported in this table 
mainly differ for the modelling of the iron feature. In the first four models the X-ray continuum is fitted with const*phabs*CompTT, 
to which a Gaussian line (Model 1), or a diskline (Model 2), or a relativistic line (Model 3) or two edges (Model 4) are used 
to fit residuals in the iron K-shell range. Model 5 includes a different Comptonization continuum model and a self-consistent, 
relativistically smeared, reflection component: const*phabs*(nthComp + rdblur * reflionx). 







Model 1: 


Model 2: 


Model 3: 


Model 4: 


Model 5: 


Component 


Parameter 


Gaussian 


Diskline 


Relline 


Two edges 


Reflection 


edge 


E edge (keV) 


_ 


_ 


_ 


7.5 ±0.1 


_ 


edge 


t(x1(T 2 ) 








6± 1 




edge 


E edge (keV) 








8.49 + »-» 

-0.07 




edge 


t(x1(T 2 ) 








6± 1 




phabs 


Nn (xW 22 cm- 2 ) 


2.2 ±0.1 


2.2 ±0.1 


2.2 ±0.1 


2.4 ±0.1 


2.7 ±0.1 


compTT 


fcr„.»,i (keV) 


0.69 ± 0.02 


0.69+™' 

-0.02 


0.68 ± 0.02 


0.62 ± 0.02 




compTT 


kT e (keV) 


3.2 ±0.1 


3.3 ±0.1 


3.3 ±0.1 


3.1 ±0.1 




compTT 


r 


14.1^* 

-0.5 


13.9 + " 

-0.3 


14.0+"- 4 

-0,2 


15.2 ±0.3 




compTT 


Norm(xl(T 2 ) 


7.6 ± 0.3 


7.4 ± 0.3 


-0.4 


8.3 ±0.3 




nthComp 


r 








1.84 + oo* 
4.9+tf 


nthComp 


kT e (keV) 










nthComp 


kT bh (keV) 










0.71 + fe 

-0.01 


nthComp 


Norm (xl(T 2 ) 










4.9 ± 0.2 


Gauss 


E (keV) 


6.57 ± 0.05 










Gauss 


a (keV) 


0.6 (frozen) 










Gauss 


Norm(xl(T 4 ) 


8.8 ± 1 










diskline 


E (keV) 




6.45 + !^ 

-0.Q7 








diskline 


Betor 




(-2.8) + ° 2 








diskline 


R,„ (GM/c 2 ) 




18i 

-6 








diskline 


R„,„ (GM/c 2 ) 




1000 (frozen) 








diskline 


i (°) 




60 (frozen) 








diskline 


Norm (xl0~ 4 ) 




9.6 ± 1 








relline 


E (keV) 






6.43 + 0.07 






relline 


Index 1 






2.8 + »- 2 

— U. 1 






relline 


i(°) 






60 (frozen) 






relline 


R m (GM/c 2 ) 












relline 


R out (GM/c 2 ) 






1000 (frozen) 






rdblur 


Betor 










-2.8 (frozen) 


rdblur 


R iB (GM/c 2 ) 










2o:f 


rdblur 


R 0Dt (GM/c 2 ) 










1000 (frozen) 


rdblur 


i(°) 










> 44 


reflion 


Fe/Solar 










1 (frozen) 


reflion 


r 










1 84 +004 


reflion 












660 (frozen) 


reflion 


Norm(xl0~ 6 ) 










9 9+0.5 
z ' z -0.4 




Total x 2 (d.o.f.) 


1489 (901) 


1463 (899) 


1464 (899) 


1519 (899) 


1463 (900) 



in addition to a th ermally compto n ized c ontinuum mod- 
ell ed with nthComp, by IZdz iarski et al. (1996), and extended 
by IZvcki et aLI dl999l) '). instead of compTT. The ^ 2 /d.o.f. was 
1555/901, without the inclusion of relativistic smearing. Since 
the iron line was found to be significantly broad in the previous 
models (Model 1, 2, and 3), we added the relativistic smearing 
using the rdblur component. The addition of this component to 
the model constitutes Model 5 in Table 3 and led to a ^ 2 /d.o.f. 
= 1466/899. The decrease of the y 1 f° r the addition of the rela- 
tivistic smearing was A^ 2 = 89 for the addition of 2 parameters 
(corresponding to an F-test probability of chance improvement 
of ~ 10~ 12 ). Assuming that iron has a solar abundance and freez- 
ing the emissivity betor index to -2.8 (value obtained with the 
diskline model) and the ionization parameter £ = Lx/(nr 2 ) to 
660 (this parameter tends to take high values but is unstable dur- 
ing the fit, that is why we preferred to freeze it), where Lx is 
the ionizing X-ray luminosity, n is the electron density in the re- 
flector, and r the distance of the reflector to the emitting central 
source, the inner radius has been estimated again to be 20 R„ and 



a lower limit to the inclination angle was found to be 44°. In Fig. 
3 we plot the residuals obtained with this self-consistent reflec- 
tion model, in comparison with the continuum and the relline 
model using the EPIC-pn data. 

The residuals found with respect to the different models de- 
scribed above do not show any evident systematic trend; the 
large x 2 could be due to mismatches in the cross-calibration be- 
tween the different instruments (for more details, see the cross- 
calibration document available on the XMM-Newton webpageQ) 
or to unresolved and unfitted features. 

4. Discussion 

We performed a spectral analysis of MXB 1728-34 observed 
by XMM-Newton on 2002 October 3rd in the 1-11 keV energy 
range. The best fit continuum model consists of an absorbed 
Comptonized component; the addition of a soft blackbody com- 
ponent does not improve the fit significantly. With respect to this 

2 http://xmm2.esac.esa.int/docs/documents/CAL-TN-0052.ps.gz 
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Energy (keV) 



Fig. 2. Top panel: EPIC-pn (black), MOS1 (red), MOS2 (green), 
RGS1 (cyan), RGS2 (blue) data points of MXB 1728-34 in the 
range 1-11 keV. Bottom panel: Residuals (Data-Model) in unit 
of sigmas for the diskline model (Model 2 in Table 3). 




3 5 10 

Energy (keV) 

Fig. 3. Top panel: EPIC-pn data points of MXB 1728-34 in the 
range 2.4-11 keV. Bottom panels: Residuals (Data-Model) in 
unit of sigmas for the continuum model reported in Table 2, for 
Model 3 including a relativistic line (relline), and for Model 5 
using a relativistic reflection component (reflion), respectively. 
Data have been rebinned for graphical purposes. 



continuum model, evident residuals are present at 6 - 8 keV, 
which can be fitted either by a relativistic line (such as diskline 
or relline) or by a self-consistent relativistic reflection model. 



4.1. The continuum emission 

The X-ray spectra of black-hole and neutron star in X-ray bina- 
ries are generally described by models that include a soft/thermal 
and a hard/Comptonized component; the electron temperature of 
the Comptonized component significantly decreases when the 
source transits from the hard to the soft state, while its optical 
depth increases. For hard spectra, Comptonization is unsaturated 
and the spectrum may be approximated by a cutoff power-law. 
In soft states, Comptonization is saturated and the spectral shape 
can be approximated with a blackbody emission at the elec- 
tron temperature. During the observation with XMM-Newton, 
MXB 1728-34 was in a low luminosity state, and presumably an 
unsaturated Comptonization is the expected spectral shape. We 
chose a continuum which consists of an absorbed comptoniza- 
tion model (comptt or nthComp) because these models gave the 
best fit to the broad-ban d (0.1 - 200 keV) BeppoSAX spe ctrum 
of MXB 1728-34 dPiraino et alj|2000l:lDi Salvo et alj|2000h . 

A different continuum was used bv lNg et aD d2010h to fit the 
sample of neutron star LMXBs. Most of the observations, in- 
cluding MXB 1728-34, were well fitted by a blackbody plus a 
disk blackbody component, absorbed by neutral interstellar mat- 
ter. The temperature of the blackbody and of the disk blackbody 
component had values between 1 .5 and 2.8 keV and between 0.6 
and 1.3 keV, respectively, except for MXB 1728-34 for which 
temperatures were found between 3 -8^ | keV and 1.9 + 0.3 keV, 
respectively. An alternative model consisting of a blackbody and 
a power-law was also used to fit the continuum. However, in 
this case the authors warn the reader about the fact that different 
continuum parameters were found when fitting the residuals at 
the Fe band with different models, meaning that line fit may not 
be realistic. Moreover, the authors note that MXB 1728-34 was 
in a low luminosity state during the XMM-Newton observation 
and that for such dim sources Compton scattering is expected 
to play an important role but was not included in their models. 
In our case, fitting the X-ray continuum with a Comptonization 
component, the parameters of the continuum do not change sig- 
nificantly when we add a gaussian or a relativistic line (such as 
diskline or relline) to the continuum. Our choice of the contin- 
uum modelling allows us to better constrain the profile of the 
broad iron line. In this context, the choice of a particular fitting 
to the line does not result in sensible changes of the parameters 
determining the continuum emission. Moreover when we tried a 
self-consistent reflection model, we found again similar values 
of the spectral parameters, which again indicates that the fit is 
stable and does not depend on the particular choice of the model 
used to fit the iron line. 

The equivalent hydrogen column inferred from the Galactic 
photoelectric absorption component, A^h ~ 2.8 x 10 22 cm -2 , 
is in agreement with t ypical values for this source (Nh ~ 
2.6- 2.7 x 10 22 c nr 2 , dDi Salvo et alJl2^|plrdnoetd] [2000: 
ID' Ai et all l2006l) V The Comptonized component can be pro- 
duced by inverse Compton scattering from relatively hot elec- 
trons {kT e ~ 3.3 keV using the compTT model or kT e ~ 4.8 keV 
using the nthComp model) off soft photons (fcr see d ~ 0.7 keV 
using the compTT model or the nthComp model). The seed pho- 
tons for Comptonization are compatible with coming from the 
neutron star surface. In fact we can estimate the size of the 
emitting region of t he soft photons, using the formula given by 
lin't Zand et al. (1999). For this, one assumes that the bolomet- 
ric luminosity of the soft photons is equal to the correspond- 
ing blackbody luminosity at the Wien temperature. The rela- 
tive gain y = 5.01 (for a spherical geometry) takes into account 
the energy gained by the photons scattered off relativistic elec- 
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trons through the inverse Comptonization process. This leads to 
a value R see d = 4.9c/s.i km, where ds.i is the distance in units of 
5.1 kpc, considering the unabsorbed flux that we extrapolated in 
the energy range 0.1-150 keV, F\, \ ~ 1.5(2) x 10~ 9 ergs cm" 2 
s" 1 . 

Usually a soft blackbody component is required to fit the 
broad band X-ray spectra of LMXBs, most frequently inter- 
preted as emitted by an accretion disc. This component is not 
significantly detected in the XMM-Newton spectrum. This may 
be ascribed to the relatively low X-ray luminosity of the source 
during the XMM-Newton observation, specially in the soft band. 
Also the bolometric X-ray flux in the range 0.1-150 keV im- 
plies a bolometric X-ray luminosity of Lx ~ 5 x lO 36 ^ 2 ^ erg 
s _1 , corresponding to 2% of the Eddington luminosity, that is 
Lpm _ 2.5 x 10 38 erg s" 1 for a 1 .4 M Q neutron star (e.g., 
van Parad iis & McClintock|[T994l) . In fact, during the high/soft 
state the disk is expected to be very close to the compact object, 
while in the low/hard state the disc should be truncated far from 
the compact object, and therefore its contribution is expected 
to be less important. We therefore conclude that the blackbody 
component is just too weak to be detected. This is in agreement 
with the results obtained from the fit of the iron feature with 
a reflection model. This indicates that the inner accretion disk 
is probably truncated far from the neutron star (R{ n > 25 km), 
and with a relatively high value for the system inclination with 
respect to the line of sight estimated at 44° < i < 60° which 
would further reduce the disk luminosity with respect to the 
Comptonized component in the hypothesis that the last one has 
a spherical geometry around the compact object. 

4.2. The iron line emission 

Recently iNg et all (1201 Oh and ICackett et aT] d2010h presented 
a spectral analysis of a sample of neutron star LMXBs ob- 
served by Suzaku and XMM-Newton, respectively, with partic- 
ular interest on the ir on discrete features in these sources . While 
ICackett et al. (2010) conclude that Fe K line profiles are well 
fitted by a relativistic line model for a Schwarzschild metric 
in most cases a nd imply a nar row range of inner disk radii 
(6-15 GM/c 2 \ lNg~et all (1201 0) conclude there is no evidence 
for asymmetric (relativistic) line profiles in the XMM-Newton 
data, although the line profiles (fitted with a simple Gaussian or 
a laor model) again appear quite broad, with Gaussian sigmas 
ranging between 0. 17 up to 1.15 keV. Another aspect relevant in 
this context is the impact of photon pile-u p on relativis t ic dis k 
lines and continuum spectra. As shown in Mill er et all (l2010h . 
while severe photon pile-up distorts relativistic disk lines and 
the disk continuum is opposing ways, the line shape does not 
sensibly depend on the pile-up fraction when its effect is mod- 
est. 

We tried different models to fit the iron line profile. Although 
a Gaussian line provides an acceptable fit of the line profile, we 
tried to physically interpret its large width using models for a 
relativistically smeared line in an accretion disk. 

Using the diskline profile (Model 2) we find the line cen- 
troid energy at 6.45 keV, compatible with a fluorescent Ka tran- 
sition from mildly ionized iron (Fe I-XX). The inner radius is 
in the range 12-21 R g . The line profile appears therefore to be 
significantly broad and compatible with a diskline profile. The 
results obtained by using a relativistic line profile correspond- 
ing to Model 3, which uses the more recent relline, instead of 
diskline, are perfectly consistent with the diskline model. The 
inner radius is estimated to be in the range 15-22 R g . 



The line profile can be equally well fitted using a self- 
consistent relativistic reflection model (Model 5). The addition 
of the rdblur component improves significantly the y 2 - This in- 
dicates that the line is indeed broad, and that the width of the 
line is in agreement with a relativistic smearing in the disk. The 
value of the inner radius is again consistent with that found us- 
ing a diskline or a relline profile, even if the uncertainty is larger 
{Ri a = 13-43 R g ). For a neutron star mass of 1.4 M Q , the inner 
disk radius is in the range 25-100 km from the neutron star cen- 
ter, and so the disk would be truncated quite far from the neutron 
star surface. The inclination angle of the system with respect to 
the line of sight is found to be > 44°, which is still compatible 
with the absence of dips in its lightcurve (which implies i < 60°). 

We also attempted to fit the iron feature using an alterna- 
tive model (Model 4), consisting of two absorption edges in- 
stead of a broad emission line (see lD'Ai et al.l l2006). Their en- 
ergies correspond to moderately ionized iron, Fe from IX to 
XVI, and highly ionized iron, Fe XXIII, respectively. The fit 
was not so good as with the previous models, since the corre- 
sponding x 1 was larger by 56 with the same number of degrees 
of f reedom. We can compare these results with those obtained 
by (iD'Aiet all [2006). using simultaneous Chandra and RXTE 
observations. The edges were found at slightly different energies 
of 7. 1 and 9.0 keV, corresponding to weakly (Fe I-V) and highly 
ionized iron (Fe XXV -XXVI), respectively. Although the fitting 
with two iron edges cannot be completely excluded yet, we think 
that the most probable explanation is that a couple of edges may 
mimic, given the relatively low statistics, the shape of a broad 
iron line. We therefore favor the interpretation of the iron feature 
as a broad and relativistic emission line produced in the accretion 
disk, since this gives a better fit of the XMM-Newton spectrum 
and very reasonable values of the reflection parameters. 
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